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1 Introduction 



The spatial string tension cr s (T) in finite temperature QCD is determined with lattice Monte 
Carlo techniques [H [2J O HI [5] by studying in asymmetric N t <C N s lattices expectation values 
of large Wilson loops with both sides in spatial directions. It is a nonperturbative quantity 
and thus not analytically calculable. However, if it is determined numerically in 3 dimensional 
(3d) SU(3) gauge theory [6], one can use the well established equations of perturbative high 
temperature dimensional reduction to compute it in full 4d finite T QCD and compare it 
with the lattice result [8j. The agreement is surprisingly good. 

Since cr s (T) is so well under control in QCD, it is a good laboratory for testing various 
models for AdS/QCD [{HEBE]- The heart of QCD computations is conformal invariance 
breaking and running of the coupling constant while the heart of well established AdS/CFT 
duality is conformal invariance of CFT. The essential part of AdS/QCD models is thus how 
conformal invariance breaking is modelled. We shall in this article apply the AdS/QCD 
model of Kiritsis and coworkers [12 \ \13 \ [T4" l I15j. approximating it in a form which permits an 
analytic solution. 



2 QCD discussion 

A first-principle method for determining the spatial string tension in hot QCD matter is to 
measure rectangular Wilson loops of size X x Y in the (x, y) plane. A potential V(L) is then 
defined by: 

V{L) = - lim hogW(L,Y) (2.1) 
and then finding a s from the large- L behavior 

L^oo Li 

Measurements are done with full 4d QCD action on a finite T lattice aNt = 1 /T -C aN s and 
they can be done for pure SU(3) Yang-Mills theory or for QCD with dynamical quarks. 

There is also a natural way to interpret the result and to derive it by splitting the problem 
in an analytic and a numerical nonperturbative part. The spatial string tension lives in the 
3d spatial space and one can also determine it for 3d SU(3) Yang-Mills theory. Since the 
only dimensionful parameter in this theory is its coupling constant gf^ of dimension mass, 
one knows that the string tension must be a number times g\ t . In fact [6l I24j. 

■2 



V^J = 0.553(1)^. (2.2) 

On a schematic level, one knows that g 2 M = g 2 (T)T, where the 4d coupling constant l/g 2 (T) 
has a simple logarithmic expansion. Thus it is natural to choose to plot 
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where the standard 2- loop expansion is inserted (with Nf = 0). The scale A CT is unknown 
but can be fitted to the data. The data for finite T SU(3) Yang-Mills theory [3] is as shown 
below in Figs. H] and [5j 
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To obtain a controlled theoretical result one at least has to compute the A a in f|2.3|) and 
not fit it to data. To do this one can start from the nonperturbative magnetic sector 3d 
number (|2.2p and work out backwards, first to the electric sector of hot QCD by restoring 
the electric mass toe ~ gT by the 2-loop relation [16] 
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and finally the hard scale ttT by the 2-loop relations [7j {(5q 
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Here g 2 (p,) is the MS running coupling in which we, for concreteness, optimized the scale 
\x so that the g 4 term in (|2.5p vanishes, i.e., // op = 47re _7E_1//22 T = 6.742 T. A thorough 
analysis of the optimisation scale dependence is given in [7j. Noting that T c = 1.15Aj^g |7J 
the logarithmic factors in l/g 2 (fi op ) could be written as log(7. 753 T/T c ), i.e., A CT in (|2.3p has 
been evaluated to be A CT = T c /7.753. 

The reliable QCD prediction then is obtained by inserting to (|2.2[) g 2 M from (|2.4p with ^ 
and m|i as given by (|2.5p and (|2.6|P 1 The prediction is shown as the continuous curve plotted 
over the range T c < T < 20 T c in Figs. H] and [5] below. The lattice data, shown in the same 
figures, extends over the range T c < T < 4.5 T c . 

One also observes that the correction terms in (|2.4p and the g e term in (|2.5p almost cancel 
each other so that to a good accuracy the result is just given by the 2-loop expression in 

([Z5D . 

The QCD computation is an entirely controlled perturbative computation, the nonpertur- 
bative part is isolated in the number (12. 2p . It works surprisingly well, perhaps even too well. 

1 Note that in a 3d gauge theory the couplings are scale independent while m 2 has both a linear 1-loop and 
a logarithmic 2-loop divergence. It so happens that the g 4 -coefHcient of the logarithmic divergence cancels for 
m% [17] and m% above is scale independent to the order shown, p,dm 2 E /dp, = 0(g c '). 
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The computation is, namely, based on integrating out the large scales ttT and gT and one 
does not expect it to work down to T c . For the pressure a similar computation shows clear 
deviations from the lattice results for T below ~ 4T C |18j . 

Clearly, to assess the reliability of the result it would be useful to further compute the 
3-loop contribution - a formidable task. 



3 AdS/QCD model 



Consider then the computation of the spatial string tension in AdS/QCD models. In bottom- 
up models conformal invariance can be broken by putting by hand some structure in the extra 
dimension z: a hard wall or some soft function of the type exp(cz 2 ) [19]. We shall use a model 
[14] the idea of which is to add a scalar field and to create this scale dynamically via the 
equations of motion. Related work on vacuum potentials is in [20] . 



3.1 Defining the AdS/QCD model 

One starts from the gravity + scalar action (in the Einstein frame and in standard notation) 
S = J^r^T {/ d 5 xV^[R~l (<V) 2 + V ((/>)] -2fd A x V^k} . (3.1) 



or, in the string frame, writing 
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One now assumes a metric ansatz 



R + 4(d^) 2 + e-l*V(^) 



ds 2 = b 2 (z) 



-f(z)dt 2 + dx 2 + 



dz 2 
/(*) 



(3.2) 
(3.3) 

(3.4) 



The four functions b(z),f(z) in the metric, the scalar field <f>(z) and the potential V{(j){z)) 
are then determined as the solutions of the three field equations following from (13. lft : 



j> 2 1 bf b 2 
6 V 2+3 b + 3 bJ = 7 v{ ^ 



4 + 3- = 0, 
/ b 

(b = b'(z), etc.) and from a fourth equation, 



(3.5) 
(3.6) 
(3.7) 



dX 



P(X) = b—, \(z) = ~ N e 



db 
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where /3(A) is the beta function of the field theory one is seeking the gravity dual for. This is 
the crucial assumption of the model. The logic here is that the energy scale of the coupling is 
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identified by 1/z, so that z — > (z — ► oo) corresponds to the UV (IR). Any similar monotonic 
function would do and actually in the definition (13. 8ft one uses E ~ b(z)/C. Eq, (|3.8j) in the 
context of a black hole horizon at some z = Zh is discussed below (see Fig. [[]). 

For given (3{X) one can from the fourth equation solve b = 6(A), then from the third 
equation / = /(A) and finally from the second equation A = X(z) so that also b(z) and f(z) 
are determined, in terms of a number of integration constants. The role of the first equation 
then simply is to fix the scalar potential. Defining 

W = -b/b 2 , b = db/dz (3.9) 

the answer simply is 

V(<f>) = 12fW 



3A 



3^W, (3.10) 



where from the /-independent Eq. f|3.6|) (which can be written in the form bW = | <fi 2 ) and 

(USD 

W(X) = W(0) exp (-1 f x dX^) , W(0) = i (3.11) 

The normalisation W(0) follows from the requirement that the boundary be asymptotically 
AdS, V(0) = 12/ C 2 . 

In spite of appearances, the potential in Eq. (|3.10p is expressible in terms of A only. A 
thorough discussion of the choice of the potential has been given in |12j-|15j. The outcome 
is that SU(3) thermodynamics is well described by the above model with the potential 

VW) = ^{ 1 + ^i^ 4/ >g(l + ^A 2 )] 1/2 } . (3.12) 

To find predictions of the model one must solve the Einstein equations numerically, as dis- 
cussed in detail on |14l 115] . For concreteness, the parameter values and the normalisation of 
A(-z) used there are 

V x = 14.3, V 3 = 170.4, X(z = 1) = 0.0242254. (3.13) 

To supplement numerical computations it is very useful to have an analytic approximation. 
Clearly in the UV A ~ g 2 is small, but the outcome of the numerical computation is that 
even at T = T c X is small, in fact ~ 0.16. Assuming A is small, A < l/y/Vs ~ 0.08, one has 

V{cf>)^^{l + V l ^V i X 7lZ ). (3.14) 



C 2 

One observes from ()3.10j) that a (3 function 

/5(A) = -A)A« (3.15) 

leads to 

V{4)) = §{ 1+ 9(F1) A?_1 + 0(A 2(9 - 1} )) • (3-16) 



2 The numerical computation includes in (|3.12|l also two more terms which enforce QCD asymptotic freedom 
in the UV. These terms are totally irrelevant for how the model in [15] describes thermodynamics of hot SU(3) 
matter. In fact, they dominate for A < 0.006 which, using pTTH]) . corresponds to T>10 58 T C 
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Thus the approximate form f)3. 15j) with 



/5o = f ^1^ = 488.8, q = f, (3.17) 

gives the approximate potential in (|3.14|) . We shall base our analysis on this approximate 
beta function. 

Analytical solutions for the metric and the scalar for the model beta function (|3. 15 j) are 
given in the Appendix. Both small z, A or large T/T c expansions and bulk thermodynamics 
can be derived from them. Physically, they are accurate enough for T>1.5T C , but clearly the 
full potential (|3.12p is required to describe the phase transition region. 

A numerical solution [151 [21] of Eqs. ([3.5p - f|3.8|) gives for z > the functions b(z), f(z), X(z) 
and the value of T = — /(z/ l )/(47r) for each value of X(z) at the horizon: A^ = A(z/ l ). Some 
properties of the functions computed are as follows ("numerics" refers to [15j): 

1. zb{z) — > C when z — > and, related to this, /(0) = 1. In numerics C = 1. 

2. There is a horizon at z = z^, f(zh) = 0- I n numerics there is a different value of Zh for 
each value of T. 

3. X(z) increases monotonically with z from A(0) = so that the leading term at small z 
is 

H*) = TTT^T = 577 > ( 3 - 18 ) 

[-(g-l)/3 log(A^)] 1/(9 ~ i:) [1140.62 log(174.67/z)] 3/7 

where A (=1/174 in numerics) is a constant fixed by the normalisation (13.13p . 

4. The entropy density of the boundary theory is given by 

° = < T ) = V 3 = 4k b3{Zh) - (3 ' 19) 

The determination of full bulk thermodynamics and, in particular, of the phase transi- 
tion temperature T c (ttT c = 1/198 = 0.88A in numerics) between the two phases (the 
high (low) T phase with / ^ 1 (/ = 1) in (|3.4p ) is discussed in [15]. 

5. The numerically computed solutions reproduce the vacuum beta function via the rela- 
tion bdX/db = /3(A) only when Zh — > oo. In fact, if there is a horizon at some z = Zh, 
the functions b(z),X(z) terminate at Zh and thus also /3(A). The outcome of a compu- 
tation of a "thermal beta function" defined by the relation /3(A) = bdX/db is shown in 
Figs. [TJand[2j One should emphasize that these are computed using the full potential 
(I3.12p . One also sees how well the approximate beta function (|3.15p reproduces the 
beta function corresponding to the potential (|3.12p . 



3.2 Deriving the spatial string tension 

The equations giving the extremal string configurations and the corresponding potential are 
well known. There are two cases: the spatial string tension at finite T is determined from a 
space-space x, y Wilson loop, the usual string tension from a time-space t, x Wilson loop at 
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Figure 1: The "thermal" (3 function, defined by the relation /3(A) = bdX/db, from numerics using 
the full form (|3 . 1 2|) of the potential. Since there is a horizon at some z = Zh, the functions b(z), X(z) 
terminate at Zh and thus also /3(A). This Zh corresponds to some T/T c , marked on the figure. The 
true "vacuum" beta function is approached when Zh — ► oo. The potential (|3 . 12[) is so constructed that 
for large A, (3 = -f A(l + 3/(8 log A)), as required by confinement. For the function at T — 1.25T C 
the asymptotic slope is — 1.1A. At T — 0.96T C the slope is — 1.47A. The QCD universal beta function 
— 6qA 2 — 6iA 3 is only approached at A = 0.006. 




Figure 2: Same as previous figure but at smaller T. For T = 0.96T C the slope is — 1.47A, already very 
close to the asymptotic large A value —3/2 • A. Even here explicitly (3{\h) > — § A/j. 

T = 0. We shall summarize the equations separately for these two cases. The end points of 
the string are always at x = ±L/2 and the string hangs in the fifth dimension so that the 
maximum depth it reaches at x = is z* < z^. We abbreviate 

b s {z) = b(z) e 2 ^)/ 3 = b{z)\ 2 '\z), b s , = b s (z*). (3.20) 

1. Space-space loop: 

The relation between L and z* is given by 

f z * dz 

L = 2 - = L(z») (3.21) 
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Figure 3: Left panel: Extremal string configurations in a conformal theory <f> = computed from 
(pOTj) for z/, = 1 and for z* = 0.6, 0.9, 0.99, 0.999. Right panel: The potential V(L) as computed from 
(|3.21[) and (|3.22|) . Three values of z* needed to reach that particular value of L are marked. 



and the potential is 



V(L) = — 
net' 



2-rra 1 



L + 



dz b s (z) 
bl 
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ira 



dz 




b 2 s (z) 
b 2 



(3.22) 



+ — / dzbj(z). (3.23) 



In (13.23D the first term separates a term proportional to L, the second term is finite when 
L — * oo and z — > (one can put e = in it) and the third term contains the singularity 
when z — ► 0; its divergence is cancelled when the energy of two independent (anti) quarks is 
included [22]. Numerical examples are given in Fig. [3] for Zy l = \. 

One observes from Fig. [3] that the large- L domain is reached for z* — > z^. To see this 
analytically, expand the integrand in (13.211) near z = z^: f(z) = f(zh)(z — Zh) + ••, b s (z) = 
b s (z*) + b s (z*)(z — z*) + .. with = z h {l — e). Then 



L 



b s (zh) 



log- 



f{zh)b s {z h ) 1 - z*/zh 



(3.24) 



Thus L diverges when — > z^ provided that b s (z} l ) < (always f{zh) < at T > 0). Since 
6 S = 6 S | A/A • [1 + 3A/(2/3)] < the inequality demands (3(\h) > — § A/j. This is true as is 
seen from Figs. [T] and [21 even down to T = 0.96T C . From (|3,23p with z* — » z^ the final result 
is 

a s = _1 7 62 ( )A 4/3 ( j_ (325) 

This is the component of the metric (|3.4p in the string frame (|3.2p multiplied by the 
string tension l/(2ita') [22]. 
2. Time-space loop. 

In this case the relations are (/* = f(z*)) 



L 



dz 



o ^/f(z)[bi(z)f{z)/(bth)-l] 



L(z x 



(3.26) 
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and 



V(L) = — 
na 



dz b 2 s (z) 



vi - raw 



blVT* L+ ^ f* dz J U lbi(z)f(z) 



2na' 



na 



(3.27) 
(3.28) 



where the z — > divergence in V can be regulated as for an cc, y loop. If / ^ 1 L does not 
diverge, the equilibrium state is that of two separate quarks [231 [19]. However, taking / = 1 
and expanding near z = z* one finds that L = (^*& s */[— ^(z*)]) 1 / 2 so that L will diverge for 
that value of z* = z m \ n satisfying 

i> s (*min) = 0. (3.29) 



From (|3.28p the string tension then is given by 



a 



2 i b (^min)^ ^ (^min)- 



(3.30) 



The functions b(z), \{z) in (|3.25|) and (|3.30p are different, in the former they are computed 
for f(z) / 1, f(z h ) = 0, in the latter f(z) = 1, z h cxd. 

3.3 Evaluating the result for spatial string tension 

The lattice results for a s in finite T SU(3) Yang- Mills we want to compare with [3] are 
measured for T c < T < 4.5 T c . The QCD computation in Section [2] should be a good estimate 
of the data even at higher temperatures, the better the larger T is, see Fig. [5j 

It is rather automatic to evaluate the expression ()3.25[) numerically [21], but to control 
various effects contributing to it we want to see how far we can get with the approximate 
potential (|3. 14[) and the approximate beta function /3(A) = — (3o\ q in f)3. 15|) . The overall 
normalisation will be fixed from the T = string tension a ~ (440 MeV) 2 , measured as 
T c / ' \[& [23] 5 using (|3,30p . This has the important property that the arbitrary normalisation 
of A cancels. 

The result f)3.25j) can be analytically approximated as follows: 
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(3.31) 
(3.32) 

(3.33) 
(3.34) 



where we first used the small-z expansions in the Appendix on the leading log level, writing in 
the arguments of the logarithms b(zh)/bo = b(zh) / (CA) = l/(Az/ l ) = nT/A, then introduced 
the explicit values q = 10/3, (5q = 489 used in the numerics [15] and finally and most subtly, 
used 0.88 A = nT c from the numerics. 
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To get the normalisation we shall use the lattice data for T c /y/d = 0.597+0.45/iV c 2 [M] and 
its prediction (13.30P in this model. To have an analytic approximation for b s (z m i n ), b s {z ra \ n ) = 
we again use the beta function /3(A) = — /?oA 9 , although the accuracy deteriorates with 
increasing A. Since b s = implies /3(A) + | A = one immediately obtains X(z m i n ) = 
(3/(2/3 )) 1/( ^ 1) and then from Eq.flHU) b(z min ) = b exp[2/(3(q - 1)]. Thus, with b = CA, 

f 3g N 2/(3(9-1)) 

b s (z m - m ) = £A(—j . (3.35) 

Numerical computation (£ = 1, A = 1/174) with the unapproximated potential (|3.12p gives 
b s {z m in) = 0.00113; (ET351) with parameter values (IXTTD gives 0.00146. From (EPUD then 

, T c A /SeX 2 ^ 1 " , , 

VW = vr£— = — — . (3.36) 



y/a ttT c \2f3 J 

With the lattice result T c /^/a = 0.6, the numerical result ttT c = 0.88A, the value of j3o and 
the correct numerical value of b s (z m i n ) (reduction of the RHS of (|3.36p by 0.00113/0.00146) 
this can be converted to 

(3.37) 



C C 6.0 
This result depends on the arbitrary normalisation of A. 

Combining (|3.32p and (|3.36[) (3q cancels and the final result is 



W = § W I < h-^ [|e(?-l)logf ] W- S > (3.38) 

Ja s tt1 c 1 - log — 



T c 1.14 



V* 1 - m lQ g" 



2/7 



49 ^ 6 T c 

A comparison of the model with lattice data is shown in Figs. [HandEl We comment on 
the result as follows: 

1. The analytic leading log approximation (I3.32P agrees excellently with the numerical 
computation using the full potential (|3.12p for T>1.6T C provided that one supplements 
it with the numerical result 0.88 A = ttT c . 

2. The parameter (3$ = ^g-Vw/Vz cancelled in (|3.39p . but if the full potential (|3.12p is 
used, there may be some dependence on it for T c < T<1.6T C . 

3. The QCD result (|2.3p was given as a numerically computed T- independent 3d quantity 
times a perturbatively computed T dependence. In analogy, here the result is given as 
the numerically computed 4d quantity T c /\fa times a T dependence computed from a 
gauge-gravity duality model. 

4. Within the range where lattice data exists, it is fitted well by both by the QCD result 
and the AdS/QCD model discussed here. For larger T the QCD prediction varies more 
rapidly with T, ~ log T/T c , than the AdS/QCD model, ~ log 2 / 7 (T/T c ). The exponent 
here is 2/(3(q — 1)) with q = 10/3, as followed from the small-A expansion of the 
potential A3. 12f) . In the extreme UV one meets with the QCD beta function with q = 2, 
corresponding to the result behaving ~ log 2 / 3 (T/T c ). If the QCD-like behavior persists, 
the present AdS/QCD model will require some modification for T>5T C . 
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4 Conclusions 



We have in this article considered a particularly simple quantity in finite T quarkless QCD, 
the spatial string tension, both in QCD and in a model for gauge/gravity duality, AdS/QCD. 
The QCD result is perfectly under control: after one nonperturbative number is determined 
numerically, the rest follows analytically via symbolic perturbative computations, which may 
be technically very demanding. For the AdS/QCD model the result needs solutions of Ein- 
stein's gravity equations and a particular contribution of this article has been developing 
analytically tractable approximations to the model in |14j . 




Figure 4: T j ^fa~ s for the range of T/T c measured on the lattice [3]. The thick points are the lattice 
data, the dashed curve is the QCD prediction (|2.3[) . the smaller points are the result of numerical 
integration of the gravity equations of this model and the continuous curve is the leading log approx- 
imation (|3.39p thereof. The last can meaningfully be extended down to T — 1.6T C . 




Figure 5: As the previous figure (which blows up the left part of this figure) but for values of T/T c 
up to 50. Note how the QCD result (dashed line) grows faster with T. 
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There are some interesting parallels and differences in the QCD and AdS/QCD results. 
Both start from a T independent nonperturbative quantity, for QCD the 3 dimensional 
Yang-Mills string tension, for AdS/QCD from a 4 dimensional computation of the T = 
string tension, expressed as T c jyfa. The T dependence in QCD comes from perturbatively 
reintroducing the scales gT and ttT, in AdS/QCD from the extra dimensional coordinate 
dependence of the metric and dilaton in the model used to break conformal invariance. 

In the range where lattice data exists, Fig. [H both calculations agree with the data within 
reasonable estimates of errors. At higher T, Fig. the AdS/QCD result, with q = 10/3 
used here, increases more slowly than the QCD result and it should be possible to distinguish 
between the two. 

There are puzzling features in both approaches. The QCD approach adds perturbative 
corrections to a nonperturbative number and it is unexpected that it works so well down 
to T c . The AdS/QCD model fits very well and elegantly bulk thermodynamics [15] using 
the potential (|3.12p . but the power of logT is quite different from that in the QCD running 
coupling (12, 7j) . This is not surprising near T c but one would expect that at least at some high 
T the matter would probe distances small enough to see QCD asymptotic freedom. 

The spatial string tension has also been measured for Nf = 2 + 1 QCD |5j. An interesting 
further project would be to study this in AdS/QCD by developing the AdS/QCD model used 
here to also include fundamental quarks. 



Acknowledgements. We thank E. Kiritsis, Mikko Laine, York Schroder and Patta Yogen- 
dran for discussions and correspondence. We are particularly indebted to Francesco Nitti 
for the Mathematica code used to solve the Einstein equations of this model and for advice 
concerning its use. 

A Appendix 

We summarize here the solution of Eqs. (|3.5|) - (|3.8|) for the beta function (|3.15p . /3(A) = — f3o\ q . 
First, from Eq.(ET8l). 

b (q - 1)/3 A9 1 

where bo is a constant, the analogue of Aqcd, the scale at which A diverges. From (|3.1ip 

C 9{q-iyiog{b/b ) 
The second equation can then be written in the form 

f = TT = -/W(W(A), (A.3) 
az az b 

from which by integration 

r 00 i / 4 \ n 
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where 



r°° 1 

T[l-n,Q]= dte'H- 71 = e~ Q — 

1 J Jq Q n 

Correct dimensions are here given by 



n n{n + 1) 
1 " Q + Q 2 



A 



(A.5) 
(A.6) 



Using (|A.3[) one can replace z as a fifth coordinate by b, A or Q. For example, 



74 / %exp 
%Jq 



-4y 



16 



9(g-l) 2 y 



T(l -n,4Q). 



9(9 -l) 2 

If the potential satisfies the constraint (|3.10p . there is a horizon with the temperature 



(A.7) 



1 

7rT 



£ 



,3Q 



dyexp[-y-16/(9(q-l) 2 y)] 



y 1 

o 



16 



bi^ n\ V 9(? - l) 2 



T(l-n,4Q(z h )) b\z h ). 



(Ai 



Using (|A.5[) one can derive various small- z (z -C 1/A, small A, large Q) approximations. One 
has, for example, 



&(*) 
1 

A^ 1 
1 

7rT 



£ 



1 + 



, 4(2 + 9(g-l) 2 ) 



9(g-l) 2 log(Az) 81(?- l) 4 log 2 (A 



(9 -l)fo log — 
Az 



4(2 + 7(g -l) 2 ) | 



1 + 



3(9-l) 2 log 2 (Az h ) 



9(9 - l) 2 log 2 (Az) 81(9- l) 4 log 3 (Az 



£ 



b(z h ) 



9(q-l) 2 log[b(z h )/b 



+ 



(A.9) 
,(A.10) 
(A.ll) 
(A.12) 



Finally, as in [14|, [T5] one can obtain the entire bulk thermodynamics by integrating s(T) 
p>(T), S (T) = b 3 (z h )/(4G 5 ). 
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